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dependence of the abelian projection reflects short distance fluctuations, and 
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I. INTRODUCTION 



Considerable progress has been achieved recently in lattice QCD as a result of combining 
two relatively old ideas: (I) improving the continuum limit behavior of the lattice action by 
adding terms that cancel the leading finite lattice spacing corrections ("Symanzik improve- 
ment" and (II) identifying a renormalized coupling which connects lattice perturbation 
theory and MC simulations |2]||. The key observation is that the disagreement between MC 
results and lattice perturbation theory can, to a large extent, be attributed to scale-independent 
(tadpole) renormalizations of the bare coupling. By converting lattice perturbation expansions 
in the bare coupling to ones using a renormalized coupling that effectively takes these tadpole 
corrections into account, lattice perturbation theory with the Wilson plaquette action can be 
reconciled with results from simulations in the scaling region J3|. Moreover, tests of asymptotic 
scaling of physical quantities are much more successful when the perturbative beta function is 
computed using such a renormalized coupling |4]-|6]J| . At the same time, the above observation 
suggests a mean-field type modification of the relation between lattice links and continuum 
gauge fields which implies that the leading order coefficients of the correction terms to 
effect Symanzik improvement have been significantly underestimated, the more so, the coarser 
the lattice f7j. Thus, besides improving lattice perturbation theory, the effectiveness of lattice 
gluodynamics simulations per se can be dramatically enhanced by working at coarse lattices 
and using a tadpole improved version (henceforth referred to as "Cornell" action 0) of the 
continuum limit improved action of Luscher and Weisz ||. The effectiveness of the method 
can be demonstrated, e.g., by computing the off-axis interquark potential at small (6 4 ), coarse 
lattices (lattice spacing a ~ 0.4 fm) in SU(3) pure gauge theory. The violation of rotational 
symmetry inherent in the Wilson action, manifest in the 0(< 40%) deviation of the off-axis 
potential from a linear-plus-Coulomb fit to the on-axis values, is reduced to 0(< 15%) using 
the continuum limit improved Luscher- Weisz action, and is essentially eliminated (of order a 
few %) by using the continuum limit plus tadpole improved action 0j. 

In the present work we apply these ideas in SU(2) pure gauge theory, focusing on con- 
finement related aspects in the framework of abelian projection (AP) picture PJTOf. After 
partial gauge fixing the original SU(N) gauge symmetry is reduced to the U(l) Ar_1 largest 
abelian (Cartan) subgroup. Under this residual group, diagonal gluon components transform 
as abelian "photons", off-diagonal gluons and quarks as doubly- and singly- charged matter 
fields, respectively. The effective abelian theory (APQCD) that results from the integration 
over the off-diagonal gluons contains a complicated assortment of abelian Wilson loop operators 
of various sizes and charges, which describe the dynamics of the abelian photons |IT|,|I"2|| , and, 
furthermore, mass terms for monopoles of different sizes and shapes These monopoles are 



identified as singularities in the gauge-fixing condition. The conjecture is then that condensa- 
tion of these abelian monopoles leads to confinement, in the spirit of the dual superconductor 
confinement mechanism in compact QED |]T^,|TB . 



One question that we wish to address in this work is whether the improvement of the non- 
abelian action leads as well to improved continuum limit behavior of the effective abelian theory 
resulting from the projection. We do this by computing the on- and off-axis potential from 
abelian Wilson loops in APQCD and comparing the violation of rotational symmetry between 
the APQCD resulting from using the Wilson action (APQCD-W) and APQCD resulting from 
using the tadpole improved action (APQCD-I). We find that the off-axis abelian potential 
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shows restoration of rotational invariance as well, allowing (at least in principle) to study the 
abelian projection in small, coarse lattices. 

The abelian projection is gauge- dependent. Early studies [|10| using local projections, e.g. 
diagonalizing an adjoint operator, did not seem to support 't Hooft's conjecture. One evidently 
successful projection is the maximal abelian (MA) [T^|T^|, corresponding in the continuum to 
D^Aff = d^A 1 ^ — igolA™, A c £] = 0, where gauge field A is decomposed in neutral(n) and 
charged(c/i) components, A^ = A n + A c ^. One nice property of MA projection is that the 



abelian monopole density is consistent with asymptotic scaling behavior in both three |T8|JT9[| 
and four [pQ| |22[| dimensional SU(2), which suggests it may be a physical quantity. However, the 
evidence in four dimensions is not indisputable, essentially because of the lack of a dimensionful 
parameter in pure gauge QCD 4 . Scaling behavior has not been observed in other projections. 

Thus, the second issue that we wish to address is whether renormalized perturbation the- 
ory and tadpole improvement allow to make more conclusive tests of asymptotic scaling of 
the monopole density. We first use renormalized perturbation theory to reanalyze the results 
which have been obtained with the Wilson action. In agreement with SU(3) results || and 
also with earlier SU(2) analyses [|,|| using the "energy" scheme coupling, the asymptotic scal- 
ing behavior of the string tension shows remarkable improvement when using the "potential" 
scheme renormalized coupling ||. In the case of the monopole density, however, although we 
do find that the degree of scaling violation is reduced when the renormalized coupling is used, 
asymptotic scaling is clearly violated for lattice spacings a > 0.1 fm, independently of the 
coupling (bare or renormalized) used. We also find evidence for (non asymptotic) scaling of 
the monopole density against the string tension, which clearly breaks down for a > 0.1 fm, as 
well. 

We then study the string tension of the improved theory (QCD-I) and the monopole density 
in its abelian projection (APQCD-I). We find deviations from scaling to be typically of the 
same order of magnitude as with the Wilson action. Monopole properties are also similar with 
those in the projected Wilson theory. In particular, we see a small scaling window for the 
monopole density in MA gauge setting in at (3 — 3.3 where the lattice spacing is quite small 
a ~ 0.1 fm, and that only for relatively large lattices (L > 12). We also find good evidence 
for (non-asymptotic) scaling of the monopole density, although not at coarse lattices. Thus, 
the improvement programme reveals that the monopole definition is plagued by certain lattice 
artifacts that do not allow to work at coarse lattices. 

A final objective of this work is to shed some light on the issue of gauge dependence of 



the abelian projection, along the spirit of previous work in three dimensions. There, it 
was shown that the difference between MA and local projections can be attributed to highly- 
correlated short distance fluctuations. Confinement, being a large-scale phenomenon, should be 
oblivious to such fluctuations. It was shown that by considering monopoles defined non-locally 
on the lattice (extended monopoles) such fluctuations are averaged over, resulting in a con- 
verging behavior between MA and local gauges. In the last part of this work, we apply similar 
considerations in the four dimensional theory, using the tadpole improved action. We find a 
similar converging behavior of the density of extended monopoles between MA and F12 gauge 
as a function of the lattice size m of the extended cube used to define the extended monopoles. 
The weakening of gauge- dependence at large physical scales is also demonstrated by showing 
that the density of extended monopoles in physical units forms a universal (independent of 
the lattice size of the extended monopole and the projection) trajectory as a function of their 
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size in physical units. These results allow some optimism that the large-scale dynamics of the 
abelian projection may after all be independent of the specific gauge used to implement it. 

The structure of this article is as follows: in Sec. II the action and the observables considered 
in this work are described. Results are presented in Sec. Ill and our conclusions appear in Sec. 
IV. 



II. METHOD 

The action used in this work is a tadpole improved version of the tree-level continuum limit 
improved SU(2) action of Liischer and Weisz. We begin by briefly summarizing the action 
improvement program. The Wilson action for SU(N) Yang Mills reads S[U] = (3 J2 p i S P i, where 

S pl = ^Re Tr(l - U pl ) = ^Tr(F? 2 ) + 0(a 6 ) . (1) 

Here = 47rao is the bare lattice coupling constant, a the lattice spacing, and we have taken 
for simplicity the plaquette U p i to be in the (1,2) plane. The continuum action is recovered by 
identifying (3 = 2N/g%. To improve the continuum limit behavior of the theory ("Symanzik- 
improvement" one adds operators that correct for the 0(a 6 ) terms. Among other possible 
choices |23f one can use 1x2 rectangular Wilson loop (labeled "rt" ) and lxlxl parallelogram 



Wilson loop (labeled u pg") terms [§] 

S[U] = (3 Cpi^Spi + (3 c r t^2S r t + (3 c pg ^S pg , (2) 

pi rt pg 

where the sums extend over all lattice points and relevant orientations of the operators. To first 
nontrivial order in perturbation theory, q = + AnaoAi, the action in Eq. @) reproduces the 
continuum action up to and including 0(a 6 ) terms, provided c pl = 5/3, c° rt = —1/12, c pg = 
(at tree level the coefficients are independent of the specific gauge group and the space-time 
dimensionality p4"|). One loop corrections A« have been computed by Liischer and Weisz for 
both SU(2) and SU(3) (Table 1 in Ref §). 

Following the convention of Ref. |7J we set q/3 — > (3i and redefine (3 = (3 p i, which makes the 
tree-level coefficient of the 1x2 term —1/20 f2"5f . Given (3 (which implicitly determines the 
strong coupling) the other two couplings are perturbatively renormalized 



A 



rt 



I - (^A pi + 12A rt ) l-n, 



20 
3 

(3 pg = -f3A pg 4ira . (3) 



As described in Ref. 0, the continuum limit behavior of the Liischer and Weisz action can 
be further improved by making the lattice links more "continuum like" . At the mean field level 
this entails setting — > Uq where one possible choice for the mean field factor uq is using 
the expectation value of the average plaquette 

uoH^ReTrf^) 1 / 4 . (4) 
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The average plaquette with Wilson action has been calculated in lattice perturbation theory 
to 0(a 2 ) ]26| and recently to 0(a 3 ) ^7|. It has also been calculated using the Liischer and 



Weisz action, Eq. (0), by Weisz and Wohlert [p4" |. However, in the latter case the result in 
numerically known to first order only: 

- log(^Re TrU pl ) = & a s , (5) 
N 2 -l f 1.72597, for N=2 

£ N = 0.366262 7T 



N \ 3.06839, for N=3 

The Liischer and Weisz action can now be tadpole improved by explicitly pulling a Uq 1 factor 
out of each link and replacing ao with a nonperturbatively renormalized coupling a s defined 
through Eq. fl5|). Since U v \ involves 4 links and U rt , U pg 6 links, one further redefines (3 — > (3uq , 
and to recover the correct continuum limit, the relative weight of the correction terms is 
readjusted by Uq = [1 — £jva s /2], using Eq. fl5|): 



p 



rt 20u 2 



1 - (jA pi + 12A rt + |0 4vra s 



3 

P pg = —^pA pg 4:iTa s . (6) 

Using Table 1 of Ref. |§ we recover for SU(3) the improved action of Ref. 0, while for SU(2) 
we find 

S = PY,S P i-Jl 1[ [l + 0.2227a s ] ^ S rt - 0.022244"* E S P9 ■ (7) 

pi Z{JU rt U pg 

The success of tadpole improvement can be seen in the value of the one-loop correction to the 
coefficient of the rectangular term, which for SU(2) becomes (1.08573 — 0.86298) — > 0.2227, a 
quarter of the original value. This is similar to what happens in SU(3), where it was shown 
that the difference between results obtained using the one-loop corrected and tree-level tad- 
pole improved actions is insignificant [0. The results reported here are in fact obtained from 
simulations using tree-level improvement only, i.e., 

S = (3^2S p i - -^-2 E Srt ■ (8) 

pi Z{JU Q rt 

For recovering the correct continuum limit it is important to realize that, in the convention of 
Ref. that we follow here, the relationship between the bare coupling ao and the simulation 
parameter (3 in Eq. (|) must be modified to take into account the absorption of c° pl in [3 |23| 

while an additional factor (1 + 1.01938 olq) is needed in the case of the one-loop corrected 
action of Eq. (|7p to account for the one-loop corrected c p i P5[ . Arguably, the RHS of Eq. (g) 
should be furthermore divided by Uq. Following Ref. f23|, this is not done here; instead, we 
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will denote the coupling resulting from such a division as a "tadpole-improved" or "boosted" 
coupling, a TI = a Wo" 4 , as in the case of the Wilson action. 

To simulate Eq. (^j) we thermalize using the heatbath algorithm, beginning with a few 
steps where u is kept fixed to 1. After obtaining a first estimate of the average plaquette (and 
therefore for w , cf. Eq. (|j)), we thermalize a few thousand times with uq computed every 
2 — 4 updates and then fed into the action, until Uq stabilizes within a few parts in 10~ 5 . For 
extracting the on- and off-axis SU(2) potential, we compute temporal T xC Wilson loops. Here 
C are spatial paths of three types: (a) straight-line spatial paths of up to 6 links, from which we 
extract the R — 1, 2 ... 6 on-axis potential, (b) planar spatial paths C — 1 x 2, 1 x 3, from which 
the off-axis potential at R = \/5, a/To is extracted, and (c) cubic spatial paths C = 1 x 1 x 1 
from which the R = potential is obtained. In the case of non-straight line spatial paths we 
sum over the possible combinations allowed given the edges of the spatial path so as to obtain 
the lowest energy (J = 0) state. Retaining only those paths minimally deviating from the 
diagonal, there are 2 such paths for (b) and 6 for (c). Measurements are separated by 20-100 
heatbath updates. 

We then perform the abelian projection to our nonabelian configurations (we henceforth 
restrict our attention to SU(2)). A lattice implementation of abelian projection was formulated 
in Refs. [0,0, in which several gauge- fixing conditions were also developed (following 't Hooft 
0). Local (generally nonrenormalizable) projections can be defined by the diagonalization 
of an adjoint operator 0]. Examples are diagonalization of a plaquette or a Polyakov line 
fT0|j . The maximal abelian (MA) projection JT6J corresponds in the continuum limit to the 



renormalizable differential gauge D^A^ = 0, where = (A 1 ^ ± iA 2 ^/\/2. Parametrizing the 
SU(2) links in the form [£|,[Lj 

/ cos^e^ sin <j> x ^ e ix ^ \ , . 

X ' A ^ - sin </> x>ft e"***.* cos0 XjA e~ ie ** J ' 1 J 

with e [0, 7r/2] and Xi @ [~ ^ 7r ]> MA projection amounts to making the transformed links 
U'* as close to the identity as possible 

max^cos(20 / x . A ) . (11) 

X,fl 

Under diagonal SU(2) transformations the phases 6 and 7 = x + 6 1 transform like abelian gauge 



field and charge-two matter field (in the continuum), respectively, while 4> remains invariant [12 
Eq. ( |TlT) is enforced iteratively; to speed convergence (typically by a factor of 3) we use the 
overrelaxation algorithm of Ref. |^9[ with parameter uj = 1.7. The iteration is repeated until 
the gauge transformation G x becomes sufficiently close to the identity at all sites 

max{l- \TzG x } < 5< 1, (12) 

with 5 = C(10~ 7 ) used as a stopping criterion. 

The abelian potential after the projection is obtained from singly charged abelian Wilson 
loops constructed from the phases 9 

W#? c = Re J J] ^ \ = cos ( E °i) = cos ^xc , (13) 

lieTxC I ieTxC 
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with the same choice of spatial paths C as for the SU(2) Wilson loops discussed above. 

Monopoles are identified in the abelian configurations using the algorithm of DeGrand and 
Toussaint |I5| . Firstly, reduced abelian plaquette angles 9 p i are defined 



dpi = 6, 



pi 



2ttN, 



pi > 



9 P i e (- 



-it, 



-7T 



(14) 



where N p i is identified with the number of Dirac strings passing through the plaquette (N p i G 
{0,±1,±2}). The net flux of monopole current out of the "elementary" (that is, of size l 3 ) 
cube C(n,fi), labeled by the dual lattice link (n, //), is equal to the sum of Dirac strings N p i 



passing through the oriented lxl plaquettes on the surfaces of the cube [10 



N, 



m=l 



P i 



(15) 



We also consider type-II |30| extended monopoles N m constructed as the number of elementary 
(m = 1) monopoles minus antimonopoles in a spatial cube of size m 3 . For the lattice density 
of monopoles we have adopted the definition of Ref. |2T[ 



Pi 



at 



£|iv m (n,4)| 



(16) 



i.e., the three-dimensional density of the time (/x = 4) components of the monopole currents, 
averaged over all time slices L t . 



III. RESULTS 

A. The abelian potential 

We first discuss the continuum limit behavior of the effective abelian theory after the 
projection. In Fig. [1| we show the on- and off-axis QCD potential with Wilson action (QCD- 
W) and the abelian potential resulting from its maximal abelian projection (MAQCD-W). 
The results have been obtained from 3100 configurations on a 6 4 lattice at /3 — 1.7, with 
measurements separated by 40-100 updates. A problem common to both Wilson and improved 
actions when working at coarse (a ~ 0.4 fm) lattices is the difficulty in establishing plateaus 
in the time direction for the correlators, since after T = 2 time slices (corresponding to 0.8 
fm) the S/N ratio has dropped dramatically. In this work we follow Ref. [[31] and evaluate the 
potential at T = 2 ||32|| . The on-axis potential is fitted to an ansatz V(r) = or — n/12r + c 
(dotted line). To set the scale we adopt the familiar practice ]33||| of using the physical 
string tension, a y /aN c =2,N f =o = y / Oph^s — 0.44 GeV, which suggests a lattice spacing a ~ 0.39 
fm, corresponding to the fit value a = 0.75. The large deviation of the off-axis points from 
the on-axis fit shows clearly the violation of rotational symmetry plaguing both the Wilson 
theory and its abelian projection at coarse lattices. Results for the tadpole improved action 
(QCD-I) and its maximal abelian projection (MAQCD-I) are shown in Fig. |^. They come 
from 3200 measurements on a 6 4 lattice at (3 = 2.4. A similar fit to the potential at T=2 



suggests a ~ 0.39 fm. A more careful estimation of the string tension using APE smearing |34 



shows that the lattice spacing has been overestimated by ~ 5%, and is rather close to 0.37 fm, 
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which is, nevertheless, sufficiently coarse for our purposes. The QCD results are in agreement 
with recent calculations in both SU(3) [[F| and SU(2) J35|, and show that the continuum limit 
behavior of the Cornell action is clearly improved, even with tree-level tadpole improvement 
only. 

The new feature emerging from these results is that not only QCD, but the abelian pro- 
jected theory as well, shows improved continuum limit behavior and restoration of rotational 
invariance, at least in MA projection. In a generic (other than MA) projection, APQCD con- 
tains abelian Wilson loop operators of various sizes and charges. In order to Symanzik improve 
such an action, the coefficients of these terms should be carefully rearranged and, possibly, new 
terms should be added. It is not obvious whether the addition of the rectangular term in the 
QCD action with the appropriate coefficient so as to effect tree-level Symanzik improvement 
will, after the partial gauge fixing, automatically fine-tune the coefficients of the APQCD ac- 
tion in such a way as to eliminate 0(a 6 ) corrections in both QCD and APQCD. In the case of 
MAQCD-W the effective abelian action is dominated by an abelian plaquette term [TT|,|r2|J5H 
and using the approximations in Ref. |36| one may argue that the observed restoration of rota- 
tional invariance is due to a corresponding improved abelian action dominating MAQCD-I. In 
that respect, it would be interesting to test the behavior of the off-axis abelian potential in a 
local projection (e.g., field-strength gauge, F12) . Although we have seen some evidence that 
rotational invariance is restored in F12 gauge as well, the F12 abelian Wilson loops are much 
more noisy and do not allow definitive conclusions to be drawn. Using anisotropic lattices may 
allow to clarify this point, as well as the issue of abelian dominance at small and coarse lattices, 
which we have not addressed here. 



B. scaling studies 

In this section we discuss the scaling behavior of the SU(2) string tension and the abelian 
monopole density. The string tension has dimension [length] -2 , while (from the physical inter- 
pretation of the monopole density as defined in Eq. (Jl6|)) the monopole density has dimension 

r l [ m ] 

[length] -3 . Thus, the monopole density in physical units reads p' m l = p\^ t a~ 3 = ^pA 3 , while 
the string tension cr = oiatoT 2 = jj^A 2 . For these quantities to be physical, the coefficient of 
A 3 in the monopole density that of A 2 in the string tension have to become constant (indepen- 
dent of a(0)) as the continuum limit is approached {(3 — > oo). In previous scaling studies (with 
the Wilson action) the bare coupling was employed to extract a scale A using the two-loop 
perturbative beta function 



51 



a(/3)A=f^V 21 expf-^ . (17) 

Before discussing the scaling behavior with the improved action, however, it may be in- 
structive to revisit the scaling behavior of the Wilson action results, this time using both bare 
(BPT) and renormalized (RPT) perturbation theory. 
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1. Wilson action 



In BPT, the bare coupling a = (vr/5) _1 is used to extract the scale from Eq. (|17|), 
while, in RPT, a renormalized coupling is employed, e.g., the "energy" coupling = 
'I — (□)) [@,|J, or the "tadpole improved" coupling ckti = cmqUq 4 Q, or, more effectively, 



2N 



t(7V 2 -1) 



the "potential" scheme coupling «y, defined from the nonperturbatively computed heavy-quark 
potential ||. In the latter case, instead of measuring the potential, one invokes the lattice per- 
turbation theory expansion of the heavy-quark potential [|37| 



V(q) 



N 2 -I4ira 



o 



2N q 2 
iV 2 -147ra v 



1 + a 



N 11 
¥ 2¥ 




NJ 

47T 



+ 0{al) 



2N 



(18) 



where J = —19.695 and —16.954, for N = 3 and N = 2, respectively, together with the 
analogous expansion of some other short-distance quantity, e.g., the average plaquette ||26"| , |27 



(-Re TrU pl ) 



el 



N 2 



2N 



-7T, 



cl a + c2 a\ + 0(al) 

c2 = (N 2 -l) (i^) 2 (o.0203 N 



1 

32 



in order to extract ay from the measured value of this quantity: 



log(-Re TrU pl ) 



N 2 -l q* 
2N t- «v(-) (1 + 5n «v) 



where 



UN, ,q\ NJ c2 cl 
5n = — log(^) + — + - + - 
67r 7T Air cl 2 



(19) 



(20) 



(21) 



At q* = tt one has 5n = —1.3386 and —0.8925, for N = 3 and N = 2, respectively. According 
to the procedure proposed by Lepage and Mackenzie for fixing the scale, q* = 3.41 0. Another 
nonperturbative coupling may by obtained by solving Eq. (|20|) to first order (i.e., ignoring Sn) 
and will be denoted by a s , so as to facilitate comparisons with the improved action [cf. Eq. (El)]. 

Consider first the density of elementary monopoles, whose asymptotic scaling behavior 
has been tested, using BPT, in Refs. [ 20 ^TJ] ) . The open circles in the inner graph of Fig. |3|. 



correspond to the 12 4 data in Fig. 1 of Ref. [20[|. Although the results are not incompatible with 
scaling, more solid evidence is required. The recently computed three-loop beta function for 
Wilson action [^HJ modifies Eq. ( p] ) by a factor (1 + 0.08324 A7ia ), which, however, does not 
improve the evidence for asymptotic scaling, since the f3 = 2.5 and (3 = 2.6 values come closer 
by less than 1%. Asymptotic scaling in MAQCD-W using the same data but renormalized 
perturbation theory is tested in the main graph in Fig. ^, with aA extracted from the ay 
coupling and q* = 3.41. Only results with the two-loop beta function are shown, since the 



(scheme dependent) three-loop correction is negligible in renormalized coupling schemes [38 



The values still do not lie on a plateau, however, the (3 dependence is significantly reduced: 
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between (3 = 2.5 (a ~ 0.085 fm) and 2.6 (a ~ 0.061 fm) the density drops by 23% using the 
bare coupling, but by much less (9%) using the potential coupling. Between (3 = 2.5 and 2.7 
(a ~ 0.045 fm) the corresponding numbers are 37% and 13% with a and ay, respectively. 
Reduced scaling violation is observed with other renormalized couplings as well, however, it 
is less effective compared to the potential scheme: between (3 = 2.5 and (3 = 2.7 we find 27% 
scaling violation with oti and 17% with a s (the "first-order potential" coupling). Results from 
larger lattices show smaller scaling violation: by reanalyzing the 16 4 data points in Table 2 
of Ref. [^1|] we find 9.3%, 4.5% and 4.2% scaling violation with a 0) a Ti and a s , respectively, 



between (3 = 2.5 and (3 = 2.6. 

The information Fig. reveals regarding the status of the MA monopole density as a 
possibly physical quantity should be properly assessed against similar tests for bona fide physical 
quantities, such as the string tension. Asymptotic scaling for the SU(2) string tension has been 
tested in bare PT J33|, as well as in renormalized PT, using the energy scheme coupling H£§. 
In Fig. |] we reproduce the BPT test (inner graph) and also show results for RPT in the ay 
scheme (main graph) |39| . As has been observed in Ref. |33| , asymptotic scaling is not satisfied 



in bare PT, although the (3 > 2.7 data from more recent calculations approach a plateau. 
Scaling violation is substantially reduced when renormalized PT is used, in agreement with 
Refs. where the qe coupling was employed. Specifically, from Table | follows that between 
a ~ 0.17 and ~ 0.4 fm the violation of asymptotic scaling is ~ 27% when the bare coupling is 
used, 16% with the tadpole improved coupling, 8% with either the energy coupling, ctE, or a s , 
the coupling obtained from Eq. fl2H| ) to first order, and is reduced to 6% with the ay coupling, 
obtained by solving Eq. (^) to second order, with q* = 3.41. 

Moreover, combining all these calculations we find behavior compatible with (non asymp- 
totic) scaling of the monopole density against the string tension, as is seen in Fig. |5] (see 
also |4T|). However, this scaling clearly breaks down at lattice spacings larger than some 
critical value somewhere between 0.12 and 0.09 fm. 



2. Improved action 

When testing asymptotic scaling with the improved action of Eq. (^), we again extract the 



scale A from Eq. ([P7D. Given the value of (3 used in the simulation, several choices of associated 
strong couplings a to be used with Eq. (|H]) are available, as with the Wilson action, e.g., the 
bare coupling, Eq. (|9|), or the "tadpole- improved" coupling ckti = aou^ 4 . The situation is 
somehow different with respect to the potential scheme coupling, ay, since, as remarked in 
Sec. II, the expansion leading to Eq. (|5|) has been (numerically) carried out to first order only. 
Thus, we do not have the analog of Eq. (|20"D whose solution (to second order) would give the 
corresponding ay coupling for the improved action. A solution to first order (in which case 
it does not matter which scale q* is the coupling extracted at) is available, of course, and 
corresponds to the nonperturbative coupling a s of Eq. (||). 

Let us now discuss the results obtained with the improved action. Starting with the density 
of elementary (m = 1) monopoles in APQCD-I, results from 8 4 , 12 4 and 16 4 lattices, for (3 
values between 2.7 and 3.7, are shown in Fig. [] (see also Table |TI|). The behavior of the 
monopole density is very similar to what has been observed in earlier studies using the Wilson 
action . In particular, the raw density of F12 monopoles appears to be independent of 
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(3 (and therefore of the lattice spacing), implying that in physical units it diverges like a -3 . 
On the other hand, the MA monopoles do seem to develop an exponential falloff which is a 
necessary condition for scaling behavior. One also notices the considerable volume dependence 
of the MA monopole density results, which should be expected if monopoles play some role in 
the confinement mechanism: the lattice spacing is a ~ 0.087 fm at (3 — 3.4 (Table |TJ) and 
the underestimation of the density on the 8 4 lattices compared with the 12 4 results reflects the 
inadequacy of the smaller lattice to provide for the typically 1 fm confinement scale (similarly 
at (3 — 3.5, where a ~ 0.075, between the 12 4 and 16 4 results). This underestimation is also in 
agreement with the expectation that it is the large monopole loops that are strongly correlated 



with confinement |4T]4"2]| 



Asymptotic scaling for the monopole density in MAQCD-I is tested in Fig. ^, using the 
bare coupling of Eq. ([5]), the boosted or tadpole- improved coupling a TI , and the "first-order 
potential", nonperturbative coupling a s . The results are, in every case, not incompatible with 
scaling, at level similar to that of the Wilson action results: between (3 = 3.4 and 3.6 (a ~ 0.087 
and 0.063 fm, respectively, corresponding roughly to the (3 € [2.5, 2.6] interval for the Wilson 
action discussed above) the scaling violation of the 16 4 data is 8.6%, 6.1% and 3.4% using 
q;o, ®Ti and a B , respectively. At even smaller lattice spacings, for fixed lattice size, the lattice 
volume becomes too small and this results in the expected "enveloping" curve |20| . 

Turning to the string tension, this is extract from linear + Coulomb chi-squared fits (without 
fixing the Coulomb coefficient to the Liischer value) to the time-dependent potential V(R, T) = 
\og(W[R, T— 1]/W[R, T]), using jackknife errors for V(R,T). It is known that the correction 



terms in the action spoil the hermiticity of the transfer matrix [43]. As a result, correlators 



show damped oscillatory behavior in T ||32|| . Using APE smearing ||34|| , we have nevertheless 
established plateaus in T in our fits of the string tension (Table |HI|) ; results from the 16 4 
lattices are shown in in Table [TV]. Above (3 = 3.4 the string tension from 12 4 lattices (not 
shown in Table |TV| ) is underestimated (commensurate with the deviation of the monopole 
density results between 12 4 and 16 4 lattices discussed above), since the lattice size is below 1 
fm in this case. This can be seen in Fig. |8|, where we test asymptotic scaling for the string 
tension using the improved action, and with the same choice of couplings as for the monopole 
density. Since the three-loop correction to the beta function is not universal, the Wilson action 
computation in Ref. |38[ in not applicable in this case. The results appear compatible with 
asymptotic scaling, especially when the boosted coupling, an, is used. Specifically, from Table 
[V| is seen that the violation of asymptotic scaling between (3 = 2.7 (a ~ 0.26 fm) and (3 = 3.6 
(a ~ 0.06 fm) is ~ 6% with axi, but more pronounced, ~ 12.5%, using either one of ao, 
a s and ag = (1 — (□))/].. 7259 (not shown). Apparently (and unlike the monopole density 
results discussed above) the nonperturbative coupling a s is not so successful here, compared 
to axi- However, this depends strongly on the value of the string tension at the coarser lattice 
considered, f3 = 2.7: the scaling violation between (3 = 2.9 (a ~ 0.2 fm) and (3 = 3.6 is 6% and 
4% with axi and a s , respectively. In general, asymptotic scaling tests using these two couplings 
give rather similar results, since the ratio a s /aTi slowly varies around 0.930 with either Wilson 
or improved action, as seen from Tables [I] and |V|. Notice, though, that the ratio does not 
monotonically decrease with /3, as does, e.g., for the ratio a s /a in the case of the improved 
action [cf. Table [V], or, in the case of the Wilson action, when one instead of using a s solves 
Eq. (|20"D to second order to obtain ay (last column in Table H). It is conceivable that knowledge 
of the second order terms in Eq. (|) and the improved action lattice perturbation expansion of 
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the potential (the analog of Eq. flUf)), which would allow to extract the «v coupling with the 
improved action, might lead to even smaller asymptotic scaling violation. 

As with the Wilson action, we can now combine the results to test (non asymptotic) scaling. 
Reading off ^fa ~ 3.8Ati from Fig. |8|(b) and estimating the scaling value for the monopole 
density at pM ~ 54A^j from Fig. 0(b) suggests that the dimensionless combination 
should saturate to 54 x 3.8~ 3 ~ 1. This is indeed verified in Fig. || By combining the 12 4 and 
16 4 results for the two quantities, their individual volume dependence cancels, to a large extent, 
in the dimensionless ratio. Thus, the evidence for scaling from this graph is more compelling 
than the evidence for asymptotic scaling in Fig. ^. It appears, therefore, that the elementary 
monopole density in maximal abelian gauge is indeed a physical quantity. Using yfa = 0.44 
GeV with pM ~ <x 3 / 2 implies a density of approximately 11 elementary monopoles per fm 3 , 
in good agreement with estimates using the Wilson action |f2TJ. However, although scaling is 



observed, its onset is at relatively small lattice spacings a ~ 0.1 fm, and therefore the situation 
for the monopole density is quite similar to what we found with the Wilson action (Fig. [5]). 

3. Comparison of Wilson and improved action 

A comparison of the scaling results discussed above with the two actions suggests the 
following: 

1. In both cases, renormalized couplings lead to improved asymptotic scaling behavior com- 
pared to the bare coupling^; among renormalized couplings, the nonperturbative ones (a s 
or ay, when available) are typically more successful than the tadpole-improved coupling 
q; T i. 

2. The string tension using the Wilson action shows more pronounced scaling violation 
compared to using the improved action, but only when analyzed with the bare coupling. 
When renormalized couplings are used, the results are rather similar for the two actions, 
for both observables tested here (string tension and monopole density), i.e., renormalized 
PT with the improved action does not improve renormalized PT analyses of Wilson 
action results significantly. Thus, the results using the improved action lead to the same 
interpretation as the one drawn from using a renormalized coupling with the Wilson 
action: 

• the string tension calculation suggests that the improvement program does allow 
confinement to be studied at relatively coarse lattices. 

• evidently a physical quantity, the monopole density in maximal abelian projection 
is, nevertheless, sensitive to small distance (a = 0.1 fm) physics. This sensitivity 
persists even when a continuum limit improved action is used. 



1 note that this is not true for the improved action if axi is interpreted as a bare coupling (see 
discussion below Eq. (0)). 
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In that respect, the merit of the tadpole- improvement program (in the form of, first, 
renormalized perturbation theory, and second, tadpole-improving the action itself), in- 
sofar as abelian monopoles are concerned, has been to point out that the asymptotic 
scaling behavior of the string tension can be extended to coarser lattices, while that of 
the monopole density can not, something that — due to the lack of asymptotic scaling 
with Wilson action and with the bare coupling — could not have been previously real- 
ized. To further clarify some of the above points, more accurate determinations of the 
string tension with the improved action using anisotropic lattices (following Ref. ||32|| ) and, 
possibly, a continuum limit improved version of Creutz ratios, should be performed ||44|| . 



C. Gauge dependence and extended monopoles 



It has been pointed out by several authors |]19 H2^|28[| that the large number of monopoles 



and the associated scaling violation observed in local, unitary gauges (such as F12 or Polyakov 
gauge) can, to some degree, be attributed to strong short distance fluctuations. Defining 
monopoles nonlocally (e.g., extended monopoles, introduced first by Ivanenko et al. ||3(J|] ) , 
may help average -at least partially- over such fluctuations and therefore reduce the gauge- 
dependence of the abelian projection. The main results of an analysis along these lines in the 



d = 3 theory with Wilson action by Trottier, Woloshyn and this author were the following [19 



1. monopoles in MA projection form a dilute plasma in d — 3. The monopole distribution 
is basically random and characterized by the average minimum monopole-antimonopole 
separation (r min ) which in physical units scales. The density of extended monopoles pj^' 
remains roughly the same for "sizes" ma < (r min ) (in physical units). 

2. In local gauges pf} t does not scale. The distribution is significantly narrower than what 
expected from a random distribution, indicating strong short-distance correlations. How- 
ever, unlike MA monopoles, drops rapidly with m, suggesting that by considering 
extended monopoles such correlations (which are absent in MA projection but create 
"spurious" monopoles in, e.g., F12) are being "washed out". Accordingly, the ratio of 
extended monopoles between F12 and MA projections is found to decrease as m increases. 

3. For extended monopoles with fixed (physical) size, p' m ' (in physical units) scales for MA 
gauge but not for local ones. However, the degree of scaling violation is found to decrease 
substantially as the extended monopole size increases. 

To what extent do such considerations apply to the four-dimensional theory? The situation 
appears quite similar on a first inspection: p]^' drops faster in the local projections than in MA, 
as can be seen by plotting the ratio = p^ 2 / Pma as a function of m (Fig. [K|). Empirically, 
we find that the extended monopole density in local gauges (which is, in d = 4, practically 
independent of the lattice spacing for all m) behaves like 

pK(0) = 4' (22) 

with C = 0.31(1) and 7 = 1.64(1), over a range of (3 values where the lattice spacing drops 
by a factor of 3 (Fig. |ll|); C and 7 depend very mildly on the particular local gauge fixing, 
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suggesting a purely geometrical origin of 7, which, however, we have not identified. The MA 
density can also be parametrized in the form of Eq. fl22"D , although the fits are not equally good. 
In this case the the coefficients decrease with (3, mildly for 7ma (Fig. |Tl|) and exponentially 
for C*ma (not shown). Since 7ma < 7fi2, the convergence between pj^and p^ 2 for fixed (3 
seen in Fig. |l^ readily follows from this parametrization. The more rapid convergence for the 
higher (3 values observed in Fig. TO is accounted for by 7F12 and 7ma being independent of (3 
and decreasing with (3, respectively. 

Moreover, one may show that, although, as in d = 3, p' m l = Cm" 7 a~ 3 (the density in 
physical units) does not scale, the degree of scaling violation decreases as a function of the 
monopole size ma in physical units. Unlike d = 3, where (3 is dimensionful and therefore this 
can be tested directly by comparing densities with m/ (3 = fixed, in d = 4 can only be implicitly 
deduced, from Eq. fl22"D ; indeed, a measure of scaling violation, as the continuum limit a — > 
is approached, is given by dp^/da e , with e < 0. From the above parametrization 



da' 



= C(ma)" 7 , e = 7-3 <0 . (23) 

ma=fixed 



Since < 7 < 3, the degree of scaling violation for the extended monopole density in unitary 
gauges will decrease as a function of the extended monopole size (in physical units), like in 
QCD 3 @. 

These results suggest that there is more contamination from short-distance physics in the 
local projection, as expected from the non-locality of the gauge condition in maximal abelian 
projection. In attempting to draw parallelisms with the d — 3 case, one should bear in mind 
that monopole d.o.f. are pointlike in d — 3 but form closed loops in d — 4. Thus, although a 
concept of minimum separation can still be devised at d — 4, it is probably not the correct way 
to describe the monopole distribution. Indeed, although we find that (r min ) is larger in MA 
than in F12 projection (by a factor ranging from 1.2 to 2 for the (3 values we have considered), 
when converted to physical units neither scales. That may explain why in MA projection 



does not remain constant over some range in Fig 10, as in d = 3, but instead starts to drop 



immediately with increasing m. The appropriate way to describe the monopole distribution is 



probably by categorizing the loops according to their length |4"2]] ; although not discussed here, 
the analog of the d = 3 case may be to examine how close to random the MA monopole loop 
length distribution found in Ref. is. 

In order to test the behavior of monopoles at large scales we plot in Fig. [12] a family of 
m = fixed trajectories for p^ 2 an d as functions of the extended monopole size ma in 
physical units, using the scale extracted from the «ti coupling. The linear behavior of the F12 
data points, with the m- independent, equal to —3, slope and the increasing like logm ordinates 
for a given abscissa, follows directly from the /^-independence of the coefficients in Eq. ( p2|) 
with 7 < 3: 



(pi 



log C + (3 — 7) log m — 3 log (ma) . (24) 



For each individual m = fixed trajectory, Fig. 12 essentially tests asymptotic scaling for p 



The most interesting feature is the formation of a universal, i.e., m-independent, trajectory for 
MA extended monopoles of large size in physical units. Unlike the F12 case, the dependence 
of Cma, 7ma on (3 (and therefore, implicitly, on a) does not allow a simple explanation of this 
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feature. The deviation of the individual m = fixed curves from this universal trajectory occurs 
at ma c A TI ~ 0.32 , 0.22, 0.16, 0.11, for m — 6, 4, 3, 2, respectively. Thus, the crossover points 
involve a common critical lattice spacing a c Axi — 0.055, which, from the scaling value ■sfo = 
3.8Ati = 0.44 GeV implies a c ~ 0.094 fm. Given that these results are from 12 4 lattices, and 
since for fixed m the continuum limit is to the left on the horizontal axis in Fig. O, it appears 
that this deviation from the universal trajectory is merely another manifestation of the finite 
physical volume effect occurring when the lattice size La is less than 1 fm, as described in 



and also observed in Sec. III. It is quite conceivable that towards the infinite volume limit 
this universal trajectory will extend to arbitrarily small physical sizes. We also notice that the 
difference between MA and F12 results decreases rapidly with increasing monopole size and 
at large physical sizes a projection-independent trajectory seems to form, indicating that also 
in d = 4 the abelian projection shows evidence of gauge independence when large scales in 
physical units are probed. 



IV. SUMMARY 

In this work the ideas of renormalized perturbation theory and tadpole improvement have 
been used to study confinement related aspects of lattice SU(2) gluodynamics. Observables 
have been studied that are either directly related to confinement, e.g., interquark poten- 
tial/string tension, or conjectured to be in the context of a dual superconductor picture, e.g., 
monopoles and abelian potential after the abelian projection. Two types of studies have been 
undertaken. Firstly, tests of the improvement program in QCD and its abelian projection, 
APQCD, specifically (a) whether asymptotic scaling and scaling is observed, and (b) what is 
the continuum limit behavior of APQCD at small, coarse lattices. Secondly, studies of ex- 
tended monopoles, using the tadpole improved action as a new, better tool, with the objective 
of shedding some light onto the issue of the apparent, albeit bothersome, gauge dependence of 
the abelian projection. The results can be summarized as follows: 

• At small, coarse lattices the degree of violation of rotational invariance is similar in Wilson 
QCD and its abelian projection. With the improved action, though, rotational invariance 
is restored in both QCD and the corresponding abelian theory, at least in MA projection. 

• Deviation from asymptotic scaling for the SU(2) string tension is observed at the 6% 
level between a = 0.06 and 0.26 fm. This order of scaling violation is quite similar with 
renormalized coupling analyses of Wilson action results, although considerably improved 
in comparison to the corresponding bare coupling analysis, even when the three-loop 
beta function is used. The quality of our calculation does not allow statements to be 
made about scaling at the 1% level. Using anisotropic lattices, will allow more accurate 
determination of the string tension for lattice spacings a > 0.2 fm. However, it does not 
seem very likely that, even with this technique, asymptotic scaling will be verified at the 
1% level for coarse and small lattices. The fact that continuum limit improvement is not 
very evident in the string tension calculation, is, however, not surprising, since the string 
tension is obtained from the standard on-axis potential. 

• Good evidence for scaling of the density of maximal abelian monopoles (using the SU(2) 
string tension to set the scale) is found. Some evidence for asymptotic scaling is seen as 



15 



well, although hampered by finite volume effects. It is, nevertheless, clear, that neither 
scaling nor asymptotic scaling hold for the MA monopole density at as coarse lattices as 
10% asymptotic scaling for the string tension does. 

• The gauge dependence of the density of abelian monopoles is significantly reduced when 
considering extended monopoles of large sizes in physical units. 

The scaling studies of the monopole density suggest that the monopole density in MA 
projection is a physical quantity, and thus seem to settle a hitherto not entirely clarified issue. 
Together with the reduced scaling violation in local projections at large scales, this may be 
considered as supporting evidence for the abelian projection picture of confinement. However, 
despite the encouraging results from the off-axis abelian potential, scaling of the monopole 
density breaks down at coarse lattices. It is conceivable that this necessitates improvement of 
the monopole identification algorithm besides improving the action p5 . 



ACKNOWLEDGMENTS 

The author would like to thank Richard Woloshyn, Howard Trottier and Pierre van Baal 
for useful discussions and suggestions, and also TRIUMF and the INT at the University of 
Washington for their hospitality during the summer of 1995. This work has been supported by 
Human Capital and Mobility Fellowship ERBCHBICT941430 and by the Research Council of 
Australia. 



16 



REFERENCES 



K. Symanzik, Nucl. Phys. B226, 187 (1983) . 

G. Parisi, in High Energy Physics- 1980, eds. L. Durand and L.G. Pondrom (American 
Institute of Physics, New York, 1981); G. Martinelli, G. Parisi and R. Petronzio, Phys. 
Lett. 100B, 485 (1981). 

G.P Lepage and PB. Mackenzie, Phys. Rev. D 48, 2250 (1993). 

F. Karsch and R. Petronzio, Phys. Lett. 139 B, 403 (1984). 

J. Fingberg, U. Heller and F. Karsch, Nucl. Phys. B392, 493 (1993). 

G. S. Bali, K. Schilling and C. Schlichter, Phys. Rev. D 51, 5165 (1995). 

M. Alford, W. Dimm, G.P. Lepage, G. Hockney and P.B. Mackenzie, Nucl. Phys. B (Proc. 
Suppl.) 42, 787 (1995); M. Alford, W. Dimm, G.P. Lepage, G. Hockney and P.B. Macken- 
zie, Phys. Lett. 361B, 87 (1995). 

M. Luscher and P. Weisz, Phys. Lett. 158B, 250 (1985). 

G. 't Hooft, Nucl. Phys. B190, 455 (1981). 

A.S. Kronfeld, G. Schierholz and U.-J. Wiese, Nucl. Phys. B293, 461 (1987). 
K. Yee, Phys. Lett. 347B, 367 (1995). 

M.N. Chernodub, M.I. Polikarpov and A.I. Veselov, Phys. Lett. 342B, 303 (1995). 

H. Shiba and T. Suzuki. Phys. Lett. 351B, 519 (1995). 

A. M. Polyakov, Nucl. Phys. B120, 429 (1977); T. Banks, R. Myerson and J. Kogut, Nucl. 
Phys. B129, 493 (1977). 

T.A. DeGrand and D. Toussaint, Phys. Rev. D 22, 2478 (1980). 

A.S. Kronfeld, M.L. Laursen, G. Schierholz, U.-J. Wiese, Phys. Lett. 198B, 516 (1987). 
T. Suzuki, Nucl. Phys. B (Proc. Suppl.) 30, 176 (1993); M. Polikarpov, Nucl. Phys. B 
(Proc. Suppl.) 53, 134 (1997). 

V.G. Bornyakov and R. Grygoryev, Nucl. Phys. B (Proc. Suppl.) 30, 576 (1993). 

G.I. Poulis, H.D. Trottier and R.M. Woloshyn, Phys. Rev. D 51, 2398 (1995). 

V.G. Bornyakov, E.-M. Ilgenfritz, M.L. Laursen, V.K. Mitrjushkin, M. Miiller-Preussker, 

A.J. van der Sijs, A.M. Zadorozhny, Phys. Lett. 261B, 116 (1991). 

L. Del Debbio, A. Di Giacomo M. Maggiore and S. Olejnik, Phys. Lett. 267B, 254 (1991). 

S. Hioki, S. Kitahara, S. Kiura, Y. Matsubara, O. Miyamura, S. Ohno, T. Suzuki, Phys. 

Lett. 272B, 326 (1991); Erratum-ibid. 281B, 416 (1992). 

M. Garcia Perez, J. Snippe and P. van Baal, Phys. Lett. 389B, 112 (1996). 

P. Weisz and R. Wohlert, Nucl. Phys. B236, 397 (1984). 

We thank Pierre van Baal for bringing this point to our attention. In some recent 



work [23,32] the 5/3 factor is explicitly retained rather than absorbed into (3. 

A. Di Giacomo and G.C. Rossi, Phys. Lett. 100B, 481 (1981). 

B. Alles, M. Campostrini, A. Feo and H. Panagopoulos, Phys. Lett. 324B, 433 (1994). 

S. Hioki, S. Kitahara, Y. Matsubara, O. Miyamura, S. Ohno and T. Suzuki, Phys. Lett. 
285B, 343 (1992). 

J.E. Mandula and M. Ogilvie, Phys. Lett. 248B, 156 (1990). 

T.L. Ivanenko, A.V. Pochinsky and M.L Polikarpov, Phys. Lett. 252B, 631 (1990). 

C. Morningstar and M. Peardon, Nucl. Phys. B (Proc. Suppl.) 47, 258 (1996). 

for a remedy of these problems, see C. Morningstar, Nucl. Phys. B (Proc. Suppl.) 53 914, 
(1997). 



17 



C. Michael and M. Teper, Phys. Lett. 199B, 95 (1987). 
M. Albanese et al., Phys. Lett. 192B, 163 (1987). 
H. Trottier, private communication. 
G.I. Poulis, Phys. Rev. D 54, 6974 (1996). 

A. Billoire, Phys. Lett. 92B, 343 (1980); E. Kovacs, Phys. Rev. D 25, 871 (1982). 

B. Alles, A. Feo and H. Panagopoulos, archive: [hep-lat/ 9609025 . 

We have used 0.37001(1) for the plaquette at (3 = 2.4 instead of the value in Table VII of 
Ref. §. 

G.S. Bali, V. Bornyakov, M. Muller-Preussker, K. Schilling, Phys. Rev. D 54, 2863 (1996). 
J. D. Stack, S.D. Neiman, and R.J. Wensley, Phys. Rev. D 50, 3399 (1994). 
A. Hart and M. Teper, Nucl. Phys. B (Proc. Suppl.) 53, 497 (1997). 
M. Luscher and P. Weisz, Nucl. Phys. B240, 349 (1984). 

an improved Creutz ratio, involving a suitable linear combination of W(R + m, R + n) 
Wilson loops with —2 < m, n < 1, \m — n\ < 2, was proposed in Ref. |2~| . 
[45] For example, one may attempt to improve the gauge condition [Eq. (|TT|)1 itself by adding 
a term involving, e.g., length-2 links so as to cancel the leading 0(a 2 ) error in the dis- 
cretization of the D^A± = continuum condition (G. Poulis and R. Wensley, in progress). 



18 



TABLES 



TABLE I. Asymptotic scaling test for SU(2) string tension with Wilson action. Different coupling 
schemes are used to get a(/3)A from Eq. (0). The bare coupling ao = i^P)' 1 , the tadpole improved 
coupling «ti = oo( D ) _1 j the energy coupling oe = ^(1 — and the potential coupling obtained 
from Eq. (p0|) with q* = 3.41, to first order, a s = — ^ log(D), and to second order, ay. Data from ||]|| 
and references therein. 



p 


< □ > 




A 


A T i 


A E 


As 


Ay 


a B ay 


2.3 


0.6024 


0.3690(30) 


62.4( 5) 


6.59( 5) 


22.29(18) 


2.49(2) 


1.12(1) 


0.936 1.199 


2.4 


0.6300 


0.2660(20) 


57.8( 4) 


6.44( 5) 


22.75(17) 


2.55(2) 


1.16(1) 


0.931 1.155 


2.5 


0.6522 


0.1870(10) 


52.3( 3) 


6.03( 3) 


22.07(12) 


2.47(1) 


1.14(1) 


0.929 1.127 


2.5115 


0.6544 


0.1836(13) 


52.9( 4) 


6.11( 4) 


22.41(16) 


2.51(2) 


1.15(1) 


0.929 1.124 


2.6 


0.6701 


0.1360(40) 


49.0(14) 


5.76(17) 


21.52(63) 


2.41(7) 


1.11(3) 


0.930 1.109 


2.635 


0.6757 


0.1208( 1) 


47.5( 1) 


5.62( 1) 


21.08( 2) 


2.36(1) 


1.09(1) 


0.931 1.104 


2.7 


0.6856 


0.1015(10) 


47.1 ( 5) 


5.62( 6) 


21.28(21) 


2.38(2) 


1.11(1) 


0.932 1.096 


2.74 


0.6913 


0.0911( 2) 


46.8( 1) 


5.61( 1) 


21.32( 5) 


2.38(1) 


1.11(1) 


0.932 1.092 


TABLE II. Runs for monopole density determination. After A^herm thermalization steps, iV measur 
measurements are taken, separated by iV separ updates. 


P 


lattice 


Athern 




A measur 


-^separ 


< 


□ > 


o [1] 
Pma 


2.7 


12 4 






100 


100 


0.5771(2) 


0.09157(27) 


2.9 


12 4 


2500 




500 


60 


0.6216(1) 


0.05295(10) 


3.1 


12 4 


2500 




500 


60 


0.6581(1) 


0.02548( 9) 


3.2 


12 4 






30 


20 






0.01641(36) 




16 4 


2500 




30 


50 


0.6730(1) 


0.01706(13) 


3.3 


12 4 


2500 




500 


60 


0.6862(1) 


0.01106( 7) 




16 4 


3000 




30 


50 


0.6864(1) 


0.01068(17) 


3.4 


12 4 


3000 




30 


70 


0.6980(2) 


0.00662(29) 




16 4 


2500 




30 


60 


0.6979(1) 


0.00724(13) 


3.5 


12 4 


2500 




50 


100 


0.7086(2) 


0.00364(22) 




16 4 


2500 




100 


50 


0.7086(1) 


0.00450( 8) 


3.6 


12 4 


3000 




30 


70 


0.7182(2) 


0.00164(14) 




16 4 






26 








0.00271(21) 
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TABLE III. String tension gt and the corresponding lattice spacing, as extracted from lin- 
ear-plus-Coulomb fits to V(R,T = fixed). 





T 




CTT+l 


a (fm) 


3.1 


3 


0.0985( 8) 


0.0965(13) 


0.141( 6) 


3.3 


4 


0.0483(11) 


0.0476(16) 


0.098(11) 


3.4 


5 


0.0378( 9) 


0.0377( 7) 


0.087(11) 



3.5 5 0.0278( 6) 0.0277( 7) 0.075( 8) 



TABLE IV. Runs for string tension determination. -/Vj ter iterations of APE smearing with pa- 
rameter cq are used. 



p 


lattice 


-^therm 


A^measur 


A SC p ar 


co 


N iteT 


< □ > 


a 


3.1 


16 4 


2500 


200 


40 


8 


20 


0.6581(1) 


0.0985(20) 


3.2 


16 4 


2500 


200 


40 


6 


20 


0.6732(1) 


0.0702(20) 


3.3 


16 4 


2500 


200 


20 


6 


20 


0.6862(1) 


0.0483(19) 


3.4 


16 4 


2500 


200 


20 


5 


20 


0.6980(1) 


0.0378( 5) 


3.5 


16 4 


2000 


200 


20 


5 


20 


0.7086(1) 


0.0278( 7) 


3.6 


16 4 




200 


20 


5 


20 


0.7184(1) 


0.0195(10) 



TABLE V. Asymptotic scaling for SU(2) string tension with improved action. Different coupling 
schemes are used to get a(/3)A from Eq. (|17|): «o = f^/?) -1 , «ti = «o( a ) _1 , and a s = — , from 
Eg- (1). 





lattice 


Ao 


Ati 




<Xs 
a 


a s 
a T i 


2.7 


12 4 


18.86(52) 


3.76(10) 


4.35(12) 


1.732 


0.935 


2.9 


12 4 


18.84(47) 


3.91(10) 


4.64(12) 


1.609 


0.936 


3.1 


16 4 


17.75(18) 


3.82( 4) 


4.72( 5) 


1.519 


0.932 


3.2 


16 4 


17.38(25) 


3.79( 5) 


4.76( 7) 


1.486 


0.931 


3.3 


16 4 


16.73(33) 


3.68( 7) 


4.68( 9) 


1.457 


0.931 


3.4 


16 4 


17.17(11) 


3.81( 3) 


4.89( 3) 


1.433 


0.932 


3.5 


16 4 


17.10(22) 


3.81( 5) 


4.93( 6) 


1.411 


0.933 


3.6 


16 4 


16.63(43) 


3.72(10) 


4.85(12) 


1.392 


0.934 
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FIG. 1. The on-(D) and off-axis (■) nonabelian (a) and abelian (b) potential, from 6 lattices at 
(3 = 1.7, using the Wilson action. 
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FIG. 2. As in Fig. | but using the improved action, at /3 = 2.4, where the lattice spacing is 
approximately the same. 
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FIG. 3. Asymptotic scaling test for the density of elementary monopoles in MA projection using 
the Wilson action (MAQCD-W). Data from Ref. |20| (circles) and Ref. [^] (squares). Open symbols: 
2-loop evolution with the ay coupling [Eq. (|20|)1 for the main graph, and the bare coupling ao = 
(7T/5) - 1 for the inner graph. Solid symbols: 3-loop evolution (bare coupling only). The dashed line 
corresponds to the scaling curve, p = [66 A] 3 , quoted from Ref. @. 
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FIG. 4. Asymptotic scaling test for the Wilson action SU(2) string tension. The scale aA is 
extracted as in Fig. |. Data from Refs. |J|.||. 
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FIG. 5. Scaling test for the density of elementary monopoles against the nonabelian string tension, 



using the Wilson action. Data from Refs. RpQ 33]. 
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FIG. 6. The density of elementary monopoles (in lattice units) in the abelian-projected improved 
theory (APQCD-I). Results shown for F12 projection on 8 4 lattices (O), and MA projection on 8 
(A), 12 4 (o) and 16 4 (□) lattices. 
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FIG. 7. Asymptotic scaling test for the density of elementary monopoles in maximal abelian 
projected SU(2) using the improved action (MAQCD-I). The scale aA is extracted (a) from the bare 
coupling «o [Eq. (§)], (b) from the boosted coupling an = oqUq 4 , and (c) the nonperturbative 
coupling a s [Eq. (pi)]. Results shown for 12 4 (o) and 16 4 (□) lattices. 
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FIG. 8. Asymptotic scaling test for the nonabelian string tension with improved action (QCD-I). 
The scale is extracted as in Fig. 0. Results are shown for 12 4 (o) and 16 4 (□) lattices. The dashed 
line in (b) corresponds to 3.8Axi- 
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FIG. 9. Scaling test for the density of elementary monopoles against the nonabelian string tension, 
using the tadpole improved action. Results are shown for 12 4 (o) and 16 4 (□) lattices. 
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FIG. 10. The ratio of densities of extended monopoles between F12 and MA projections using the 
improved action, as function of the size m of the lattice cube defining the type-II extended monopole. 
Results from 12 4 lattices at p = 2.7 (□), f3 = 3.1 (O), f3 = 3.3 (o), = 3.5 (A), 
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FIG. 11. Fitting the lattice density of type-II extended monopoles according to Eq. (]22): CV12 
(o), 7F12 7ma (■)• Dashed line corresponds to C=0.31 and dotted line to 7=1.64 (see text). 
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FIG. 12. The density of extended monopoles as a function of their size for F12 and MA projections. 
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